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We investigate by numerical simulations the pattern formation after an oscillating attractive- 
repulsive obstacle inserted into the flow of a Bose-Einstein condensate. For slow oscillations we 
observe a complex emission of vortex dipoles. For moderate oscillations organized lined up vortex 
dipoles are emitted. For high frequencies no dipoles are observed but only lined up dark fragments. 
The results shows that the drag force turns negative for sufficiently high frequency. We also suc- 
cessfully model the ship waves in front of the obstacle. In the limit of very fast oscillations all the 
excitations of the system tend to vanish. 
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I. INTRODUCTION 

The realization of Bose-Einstein condensate (BEC) in atomic gases have boosted intense theoretical and experimen- 
tal investigation about its exotic properties. BEC is a paradigm of a quantum fluid and in weak interacting case can 
well be described by the Gross-Pitaesvkii equation pQ . Eventually BEC spread to other systems like exciton-polaritons 
[2 , offering new possibilities for experimental tests. An interesting feature of a quantum fluid is its contrasting be- 
havior as compared with a classical fluid. The flow of a quantum fluid past an obstacle only generates drag force 
above a certain subsonic critical velocity and the energy can be dissipated into collective excitations of the fluid. This 
dissipation can be inferred from numerical experiments by the mean drag on the obstacle [3 . A superfluid behavior 
is revealed below this velocity where nucleation of vortices never occur and no excitations are generated [3H9] . For an 
appropriate velocity and size of the obstacle, a Benard-von Karmam vortex street can be generated [TP] , 

There is also a supersonic critical velocity where oblique vortex streets are transformed into stable oblique dark 
solitons [11]. For higher velocities, the general picture of the diffraction pattern in the supersonic flow past a disk- 
shaped impenetrable obstacle consists of two different parts separated by the Mach (or Cherenkov) cone [T2] . Outside 
the Mach cone there is a region of linear waves that we will refer them simply as ship waves [13] [14]. Inside the 
Mach cone a pair of oblique dark solitons is gradually formed behind the obstacle if the radius of the obstacle is 
of healing length order. For greater radius more pairs of oblique solitons are generated. Interaction of solitons was 
studied in p~5j [16] where it was found that the angle between dark solitons decreases as the obstacle radius increases 
for a fixed supersonic velocity of the flow. In previous experimental works [I3j [T71 [18] the existence of such nonlinear 
structures were suggested. However, only recently the generation of stable oblique dark solitons was experimentally 
demonstrated in the flow of a Bose-Einstein condensate of exciton-polaritons past an obstacle p~9j [20] . A numerical 
study to support these experimental findings was done in [21], and the observation of vortex dipoles in an oblate 
atomic Bose-Einstein condensate [22] suggests that the supersonic studies can also be carried in this system. 

In atomic BEC, obstacles are typically represented by detuned lasers that can be effectively be attractive (red- 
detuned) or repulsive (blue-detuned) obstacles, by the use of Feshbach resonances. The first numerical study of 
attractive obstacles was done in Ref. [23] where it was established the critical velocity to the formation of vortices and 
corrects the velocity found in [3] in the case of repulsive obstacles. Numerical studies revealed that turbulence can 
also be achieved and studied by spatial oscillation of a repulsive obstacle [24] . A clever way to control the formation of 
vortices moving attractive and repulsive laser beams was proposed in [25] . There is a special form of moving potential 
that no radiation is generated at supersonic velocities [26] . The disappearance of gray soliton and phonon excitations 
was demonstrated in [27] by oscillating a repulsive obstacle in a quasi- ID trapped BEC at high obstacle velocities. It 
was found in [28] that vibration of an obstacle modulates the vortex street, theoretically predicted in [TO] , breaking 
a symmetry. 

In the present work, we study the flow of a BEC past an oscillating attractive and repulsive obstacle. The motivation 
is to answer the question can we get rid of drag for very fast oscillations? We investigate different regimes from slow 
to very fast oscillations. Since we are working in the supersonic regime we can divide the study inside and outside 
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the Mach cone as follows. 



II. MODEL EQUATIONS 

We consider the flow of an atomic Bose-Einstein condensate (BEC) past an obstacle in the framework of Gross- 
Pitaevskii (GP) mean field approach. In the rest frame, the condensate is well described by the macroscopic wave 
function \I/ = ty(x,y,z,t) obeying the time-dependent GP equation 

iti— = V 2 ^ + U ext * + * 2 ^, 1 

ot 2m m 

where V 2 = d 2 + d 2 + <9 2 , the external potential U ext = U tr ap(%, y,z) + U(x + vt, y, z, fit) is represented by the sum of 
a harmonic trap U tra p and a time-dependent obstacle potential U that oscillates with frequency fi, m is the atomic 
mass and a is the s-wave scattering length. 

We will limit our study to the case of the quasi-2D limit, i.e., we have a strong harmonic confinement in the z 
direction. In this regime we can approximate ^(x, y, z, t) = ifj(x,y,t)(j)(z)e~' lfJ ' z / h , where <j)(z) and \i z are the ground 
state and energy respectively for the confinement in z direction [25j[29]. Substituting in Eq.Q and integrating in z 
direction we obtain 
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^ = ~L {d * + d v^ + U ^ + SMV, (2) 

where g = 47raft 2 ra -1 j (j) 4 (z)dz is the effective interaction in two-dimensions. We consider here that the obstacle 
runs close to the center of the trap. In this region the condensate is almost homogeneous and the potential in x and 
y directions is weak as compared to the obstacle potential. So the harmonic potential is neglected for studying the 
excitation caused by the obstacle. 

We introduce dimensionless variables x = x/£, y = ?//£, t = gn^t/h, ip = tf/y/rio, U = U/gno, fl = flfi/gn^, the 
Mach velocity M = v/c s , where no is a characteristic 2D density of atoms at the center of the trap, £ = h/ \frfm^g 
is the characteristic length and the sound velocity c s = H/m£. Typical experimental values are £ ~ 0.3/mi and 
h/gno ~ 0.18ms [22 . Thus for fi ~ 1 we have oscillations of the order of kHz well within experimental reach. 
Substituting in Eq. Q and after dropping the tildes for convenience we get 

i d -^ = -\(dl + dl)^ + U^+\^, (3) 



where U = U{x + M£, fit). The energy is given by 



E(t) — J dxdy 



i|v X) ,vf + WI 2 + ^l 4 



(4) 



and the rate of energy in time is given by 



Turning off the oscillation (ft = 0) the second term vanishes and we identify the first term as M times the drag force. 
Also when M = only the second term is responsible for the excitation of the system. 

For computational purposes, in Eq. ^ we make a global phase transformation ip' = e lt ip and later a Galilean 
transformation x' = x + Mt, t' = t leading to 



^ = + d 2 y W - iMd x ip - V> + + C/V , (6) 



where U = U{x, y, fit), the primes were omitted for convenience and subscripts here means derivatives. This equation 
describes the system in the obstacle reference frame. 

The obstacle is a laser beam that continuously oscillates from blue-detuned to red-detuned and vice-versa, which 
can be written as 



U (x, y : fit) = Uq cos (fit) exp 



-2 (x 2 + y 2 



(7) 
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where Uo and wq are the amplitude and the beam waist of the laser, respectively, and CL = 2tt/T is the oscillation 
frequency of the detuning in a period T. 

III. INSIDE THE MACH CONE 

For CL = Eq. (|6| supports stable oblique solitons [TTJ [I2j [16] in the obstacle frame as 

1>{*,v) = 1 + eC , (8) 

where e ta = A + iv, ( = 2u[xsm0 — ycosO], v = yl — A 2 , A = Msin#, and is the angle between the soliton 
and the horizontal axis. Msin# = ±1 defines the Mach cone and thus solitons can be found only in the region 
-arcsin(l/M) < < arcsin(l/M). 

We have solved the Eq. ([3| numerically in the supersonic regime using Uo = 25 and wq = 1. In Fig. [I] we show the 
results for supersonic flow for different oscillation frequencies. Fig. [T] depicts the case of CL = where we reproduce 
the formation of oblique dark solitons [12 j. Outside the Mach cone there is a stationary wave pattern created by 
interference of linear waves. Inside the Mach cone there are two oblique dark solitons that decay at the end points 
into vortices, situated symmetrically with respect to the direction of the flow. As we turn on the oscillations observe 
the emission of dark fragments. For CL = 0.5 these fragments can be identified as vortex dipoles and form a pattern 
of "5 in a dice". As the frequency is increased to CL = 1.5 these fragments stand well aligned as vortex dipoles as can 
be identified by the phase plot (see Fig. [2|. These dipoles are followed by a secondary radiation emission, identified 
as a straight line almost parallel to the Mach cone. As the frequency is further increased to CL = 10 the fragments 
can no longer be identified neither as a single vortex nor as vortex dipoles. By looking at phase the fragments are 
identified as short gray solitons that propagate obliquely to the flow, analogous to the ones observed in Ref. [30]. To 
check the (non)vorticity character after some time of fragments formation we turned off the intensity of the obstacle 
and no decay into vortices were observed. 
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FIG. 1: Pictures of the diffraction pattern at fixed time t — 100 with the set of parameters: Mach number M — 3, potential 
intensity Uo — 25 (the laser beam starts repulsive) and width wo = 1 for different frequencies Q = 0, 0.5, 1.5, 10 and 20. 

One can explain the general behavior as following. For Q = 1.5 the oscillation acts as a "chopper" that turns on 
and oif the dipole emission. In this specific case the on time is more than enough to generate vortex dipoles and thus 
we have excess of energy that is ejected as secondary radiation. In the fast oscillating regime (CI = 10) the time the 
oscillation is on is not enough to form dipoles and just small dark solitons can be seen. As the frequency is around 
Cl = 20 practically no more fragments can be seen. To check the consistency of our analysis we studied the number 
of fragments as a function of CL. One can estimate that rate of fragments emission is close to 1, meaning that at each 
period one fragment is emitted. The linear behavior confirms the modeling of the oscillating obstacle as a "chopper" . 

Thus we identified four regimes, for CL ~ 0.5 we have strange patterns as the "5 in a dice", for CL ~ 1.5 we have lined 
vortex-dipoles with secondary radiation, for CL ~ 10 the dipoles are suppressed and give place to small dark solitons 
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FIG. 2: Left panel: Phase of the diffraction pattern that corresponds to Fig.[T]for frequency Q = 1.5 which provides us evidence 
of vortex dipole formation. Right panel: inset showing vector velocities of the vortex dipole in the selected region. Note that 
the velocity in x direction was not shifted by the constant flow velocity, M, for visualization purposes. 



that propagate obliquely to the flow, finally for Q ~ 20 practically no excitation can be seen inside the Mach cone. If 
the laser beam starts red-detuned (attractive) we have similar results. 

Next we develop the analytical theory to describe the remaining ship waves located outside the Mach cone and 
compare it with numerical simulations. 

IV. OUTSIDE THE MACH CONE 

Ship waves are formed in front of the obstacle. The theory for a non-oscillating obstacle was previously studied for 
a S function in ref. [31] , where it was found that the density changes are given by 



with 



and 



q(k,r, M) 



Sn = Vo q(k, r, M) x cos (kr cos fi — ^ 

2fc [(M 2 - 2)k 2 + 4(M 2 - 1)]V4 



7TT [(M 2 - 2)k 2 + 6(M 2 - 1)] V2 ' 



k = 2y / M 2 cos 2 // - 1. 
The angles \i and rj are defined according to Fig. [3j and Eq. (11) is valid if 



1 \ I 

arccos ( — < ri < arccos — : 
\M J ~ 1 ~ \M 

so that the linear waves exist only in the region outside the Mach cone. 




(9) 



(10) 



(11) 



(12) 



FIG. 3: Left panel: Coordinates that defines the radius- vector r and the wave vector k. The latter one is normal to the wave 
front which is shown schematically by a curved line. Right panel: Numerically calculated wave pattern for a fast oscillating 
obstacle at fixed time t — 100 with the set of parameters: M — 2, Uq — 25, wq — 1 and Q = 40. Solid line (red) corresponds to 
linear analytical theory (Eq. 13 ) for the line of constant phase. 
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According to [31 j, one can find the shape of the lines of constant phase (wave crests) <E> = kr cos \i in a parametric 
form 



4 $ 

x = r cos x = Tip cos r/(l — M cos 2rj), 
~Y sinr/(2M 2 cos 2 77 — 1), 

rv 



(13) 



y — r sm \ 

so that for small values of 77, corresponding to waves in front of the obstacle, these lines take a parabolic form 



(M 2 - l) 3 /2 
KU) 2VM 2 - 1 (2M 2 - 1)$" " 



(14) 



Predictions of the analytical theory are compared with the numerically calculated wave pattern in Fig. [3] and excellent 
agreement is found. So, the theory previously developed in [31] remains valid even for a fast oscillating obstacle. 

For a fast oscillating we assume that the resulting ship waves can be computed by the Huygens principle, i.e., by 
the superposition of stationary densities generated by obstacles at different positions along the flow. Averaging over 
a period this can be expressed as 



1 / 7T 

5n osc = Vq q(k, r, M) x — J cos(Qt) cos ykr cos fi — — + kMt 



dt, 



(15) 



where the term kMt was added representing phase change due to the obstacle movement along the time. After 
integration in time one obtains 



8n osc = -V oq (k,r,M) ( ^3^2 ; 



j sin ^ 



kr cos fi — ^ + kMT^J — sin (kr cos fi — ^ . 



(16) 



In the region in front of the obstacle where y = (i.e. 77 = 0) and x < 0, the wavelength A = 2tt /ko is constant 



with ko = 2\jM 2 — 1. Therefore x(0) = — rcos/i and the perturbations of the condensate density take the form 

7T 2tt/c M\ 



Sn osc = -V oq (k ^M)l 2 _ 



sm 



-kox 



- 1 - sm \-k Q x- -J 



(17) 



where r = \x\. The above formulae shows that increasing Q the magnitude of the Sn osc decreases as 1/fi and in the 
limit of Q — >• 00 the ship waves vanish. The plot illustrating this behavior is shown in Fig. [4j As we see, Eq. (17) is 
accurate enough almost everywhere, for Q > 40, except in the small vicinity of the obstacle. Although the obstacle 
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FIG. 4: (a) and (b) Profile of intensity in front of the obstacle for x < 0, y = with the set of parameters: M = 2, Vq — 2, 
Uq — 25 (the laser beam starts repulsive) and wq — 1. Solid lines (black) correspond to linear analytical theory, Eq. (17), and 
dashed lines (red) to numerical solution of Eq. (|3}. (c) The solid line corresponds to the higher magnitude value of the 5n osc 
close to x — —40 and the red circles correspond to this magnitude calculated numerically at the same position. 



in the theory is represented by a delta function, our numerical simulations using a narrow Gaussian potential as the 
obstacle provide results in very good agreement with our extended theory. 



V. DRAG FORCE 
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We also computed the drag force in the x direction as 



F x (t) 



L 



(18) 



where A defines an infinite region of the fluid around the obstacle. For practical purpose we took the integration 
along our whole grid. In Fig. [5] we show the average drag taken at one period of oscillation. For slow oscillation 
frequency ft < 5 we observe that drag is decreasing and positive as expected since both non-oscillating attractive and 
repulsive potentials causes positive drag [25] . However, for Q > 5 the drag is always negative and vanishes in the 
limit of Q — >• oo. So, the answer for the question can we get rid of drag for very fast oscillations initially proposed 
is yes. Surprisingly, the mean drag also vanishes at a small region of low frequencies and this is a non-intuitive and 
remarkable result. 




FIG. 5: Upper panel: Average drag force in the x direction, < F x >, as a function of the frequency Q for different potential 
intensities. Middle panel: response function amplitude A\ as a function of ft. Lower panel: relative phase Si as a function of 
Q. Here we can see that the sign of < F x > depends on the relative phase between the forcing potential and the response main 
mode, Eq. pi). 



The system can be seen as a forced oscillator. In our case the oscillating potential forces the system and we obtain 
as output an oscillating |^| 2 . Thus the drag force in the x direction can be explicitly written as 



F x = U cos(nt)R(t), 



where R(t) is a response function given by 



R(t) 



9 / x exp 
< J A 



-2(x 2 +y 2 ) 



\ip(x, y, t)\ 2 dxdy . 



(19) 



(20) 



We observed numerically that R(t) is periodic with period T and thus can be written as a Fourier series as 

R(t) = A + A x cos(Ot + ft) + A 2 cos(2ft£ + S 2 ) H , (21) 

where A' s are amplitudes and 8' s are relative phases to the forcing potential. Averaging the drag force in time we 
have 



<F X >= 



U cos(Qt)R(t)dt , 



and only the second term of the series survives giving 



<F X >-- 



UqA 
ft 



- cos(5i) . 



(22) 



(23) 
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Thus the sign of < F x > depends on the relative phase between the forcing potential and the response main mode. 
We computed 



and from < F x > and / we obtained A\ and S\. 

A negative drag can be interpreted as a force in the upwind direction, meaning propels the laser. The question of 
the energy balance can be explained from the oscillating laser. As it attracts and repels the condensate it pumps 
energy to the system that causes an upwind force to supersede the downwind force due to the movement of the laser. 
This upwind force is only generated in the moving and oscillating obstacle. For standing (M = 0) oscillating obstacle 
the system is radially symmetric and no drag is generated. 

VI. CONCLUSIONS 

We have studied the wave pattern generated by an oscillating obstacle in the supersonic flow of a quantum fluid. 
Turning on oscillations causes disruption of the oblique solitons into dipoles. For Q = 1.5 the dipoles are emitted 
organized as a vortex dipole street. For increasing frequencies dipoles change gradually orientation in clockwise 
direction and their bunch resembles the oblique solitons. Finally for very high frequencies the angle of emission 
increases and vortices vanishes. For the waves in front of the fast oscillating obstacle, we could further extend the 
model previously developed for non-oscillating obstacle. These waves were shown to gradually diminish according to 
the averaging of emission of linear waves out of phase. Combined results, both ship waves and soliton tend to vanish 
for high frequencies leading to a vanishing drag. Remarkably, the mean drag also vanishes at a small region of low 
frequencies during his change of sign from positive to negative. So, even a very powerful laser fast oscillating from 
red to blue-detuning could pass through an atomic BEC without generating vortices or solitons. This result could be 
experimentally checked with existing setups [22]. Analogous experiments could also be performed with condensates 
of exciton-polaritons [19]. 
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